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In this talk I review the calculation of the third order corrections to the heavy quarkonium spectrum in 
the nonrelativistic effective theory framework and its application to the phenomenology of top quark threshold 
production. 



1. INTRODUCTION 

The theoretical study of the nonrelativistic 
heavy quark-antiquark system ^ and its ap- 
plication to bottomonium |^ and toponium 
physics is of special interest because it relies en- 
tirely on the first principles of QCD. The sys- 
tem in principle allows for a perturbative treat- 
ment with the nonperturbative effects being well 
under control and with no crucial model depen- 
dence. This makes the heavy quarkonium to 
be an ideal place to determine the fundamental 
QCD parameters such as the heavy quark mass 
rUq and the strong coupling constant as- An es- 
sential progress has been made in the theoretical 
investigation of the nonrelativistic heavy quark 
dynamics based on the effective theory approach 
Hj. The analytical results for the main parame- 
ters of the nonrelativistic heavy quark-antiquark 
system are now available up to next-to-next-to- 
leading order (NNLO) (sec ||] for a brief review). 
The NNLO corrections have turned out to be 
so sizeable that it appears to be indispensable 
also to gain control over the next-to-next-to-next- 
to-leading order (N'^LO) both in regard of phe- 
nomenological applications and in order to un- 
derstand the structure and the peculiarities of 
the nonrelativistic expansion. The calculation of 
the heavy-quarkonium spectrum to 0{a^) [^,0 
was the first breakthrough in the N^^LO analy- 
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sis. Recently, the ©(a^lnQ^) corrections to the 
heavy-quarkonium production and annihilation 
rates have been also obtained The subject 
of this talk is restricted mainly to the analysis of 
the spectrum. 

In the next section we briefly introduce the 
basic ingredients of the nonrelativistic effective- 
theory formalism. In Sect. 3 the structure of the 
third order corrections to the heavy quarkonium 
ground state energy is outlined and the final re- 
sult is presented. In the last section this result is 
applied for the phenomenological analysis of the 
top-antitop pair production near the threshold. 

2. EFFECTIVE THEORY OF NONREL- 
ATIVISTIC HEAVY QUARKS 

The nonrelativistic behavior of the heavy- 
quark-antiquark pair is governed by a compli- 
cated multiscale dynamics. In the nonrelativis- 
tic regime, where the heavy-quark velocity v is 
of the order of the strong-coupling constant as, 
the Coulomb effects are crucial and have to be 
taken into account to all orders in a^. This makes 
the use of the effective theory mandatory. The 
effective-theory approach allows us to separate 
the scales and to implement the expansion in v 
at the level of the Lagrangian. Let us recall that 
the dynamics of a nonrelativistic quark-antiquark 
pair is characterized by four different regions and 
the corresponding modes (i) the hard re- 

gion (the energy and three-momentum scale like 
iTLg); (ii) the soft region (the energy and three- 
momentum scale like mqv); (iii) the potential re- 
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gion (the energy scales like mqv'^, while the three- 
momentum scales like niqv); and (iv) the ultrasoft 
region (the energy and three-momentum scale like 
rugv^). Nonrelativistic QCD (NRQCD) g is ob- 
tained by integrating out the hard modes. Sub- 
sequently integrating out the soft modes and the 
potential gluons results in the effective theory of 
potential NRQCD (pNRQCD) 0, which con- 
tains potential heavy quarks and ultrasoft glu- 
ons, ghosts, and light quarks as active particles. 
The effect of the modes that have been integrated 
out is two-fold: higher-dimensional operators ap- 
pear in the effective Hamiltonian, corresponding 
to an expansion in v, and the Wilson coefficients 
of the operators in the effective Hamiltonian ac- 
quire corrections, which are series in as- 

The theory of pNRQCD is relevant for the de- 
scription of the heavy-quarkonium system. Let 
us recall its basic ingredients. In pNRQCD, 
the (self)interactions between ultrasoft particles 
are described by the standard QCD Lagrangian. 
The interactions of the ultrasoft gluons with the 
heavy-quark-antiquark pair are ordered in v by 
the multipole expansion. For the N^LO anal- 
ysis, only the leading-order (LO) emission and 
absorption of ultrasoft gluons have to be consid- 
ered. They are described by the chromoelectric 
dipole interaction. The propagation of the quark- 
antiquark pair in the colour-singlet and colour- 
octet states is described by the nonrelativistic 
Green function of the effective Schrodinger equa- 
tion. The LO approximation for the Green func- 
tion is given by the Coulomb solution, which 
sums up terms singular at threshold and describes 
the leading binding effects. The corrections due 
to higher-order terms in the effective Hamilto- 
nian can be found in Rayleigh-Schrodinger time- 
independent perturbation theory as in standard 
quantum mechanics. 

Let us now turn to the problem of perturbative 
calculations in the effective theory. Both NRQCD 
and pNRQCD have specific Feynman rules, which 
can be used for a systematic perturbative expan- 
sion. However, this is complicated because the 
expansion of the Lagrangian corresponds to a par- 
ticular subspace of the total phase space. Thus, 
in a perturbative calculation within the effective 
theory, one has to formally impose some restric- 



tions on the allowed values of the virtual mo- 
menta. Explicitly separating the phase space in- 
troduces additional scales to the problem, such as 
momentum cutoffs, and makes the approach con- 
siderably less transparent. A much more efficient 
and elegant method is based on the expansion by 
regions |], which is a systematic method to ex- 
pand Feynman diagrams in any limit of momenta 
and masses. It consists of the following steps: (i) 
consider various regions of a loop four- momentum 
k and expand, in every region, the integrand in 
Taylor series with respect to the parameters that 
are considered to be small there; (ii) integrate 
the expanded integrand over the whole integra- 
tion domain of the loop momenta; and (iii) put 
to zero any scaleless integral In step (ii), dimen- 
sional regularization, with d = 4 — 2e space-time 
dimensions, is used to handle the divergences. In 
the case of the threshold expansion in v, one has 
to deal with the four regions and their scaling 
rules listed above. In principle, the threshold ex- 
pansion has to be applied to the Feynman dia- 
grams of full QCD. However, after integrating out 
the hard modes, which corresponds to calculat- 
ing the hard-region contributions in the thresh- 
old expansion, it is possible to apply step (i) to 
the diagrams constructed from the NRQCD and 
pNRQCD Feynman rules. Equivalently, the La- 
grangian of the effective theory can be employed 
for a perturbative calculation without explicit re- 
strictions on the virtual momenta if dimensional 
regularization is used and the formal expressions 
derived from the Feynman rules of the effective 
theory are understood in the sense of the thresh- 
old expansion. In this way, one arrives at a formu- 
lation of effective theory with two crucial virtues: 
the absence of additional regulator scales and the 
automatic matching of the contributions from dif- 
ferent scales |6|pH (see also[|l^). The second 
property implies that the contributions of differ- 
ent modes, as computed in the effective theory, 
can be simply added up to get the full result. 

3. HEAVY QUARKONIUM SPECTRUM 

AT Oia^rrig) 

The analysis of the heavy quarkonium spec- 
trum at 0{a^mq) i.e. the third order correc- 
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tions to the Coulomb approximation involves two 
basic ingredients: the effective Hamiltonian and 
the retardation effect associated with the emis- 
sion and absorption of dynamical ultrasoft glu- 
ons. The general form of the Hamiltonian valid 
up to N'^LO is given in A detailed discussion 
of the effects resulting from the chromoelectric 
dipole interaction of the heavy quark-antiquark 
pair to the ultrasoft gluons relevant for pertur- 
bative bound-state calculations at N^LO can be 
found in §Q. 

For vanishing angular momentum we can write 
the energy level of the principal quantum number 



(1) 



where the leading order Coulomb energy is = 
—C'pa'^mq/'in^. SeI^^ stands for corrections of 
order a^. The first and second order corrections 
can be found in p^ , p"5|Jl^ for arbitrary n. 

The O (a^) corrections to the energy levels 
arise from several sources: (i) matrix elements of 
the N'^LO operators of the effective Hamiltonian 
between Coulomb wave functions; (ii) higher iter- 
ations of the NLO and NNLO operators of the ef- 
fective Hamiltonian in time-independent pertur- 
bation theory; (iii) matrix elements of the N^LO 
instantaneous operators generated by the emis- 
sion and absorption of ultrasoft gluons; and (iv) 
the retarded ultrasoft contribution. 

The contributions to En for vanishing j3- 
function were computed for general n in |^ . The 
corresponding correction to the ground state en- 
ergy gets contributions from all four sources and 
reads 
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where the color group invariants are Cf = 4/3, 
Ca = 3, Tf = 1/2, ni is the number of light- 
quark flavors, 5 is the spin quantum number 
and La^ = — ln(Ci?as). The numerical value of 
the QCD Bethe logarithm Lf = -81.5379 can 
be found in p|,p^. MS scheme for the renor- 
malization of asip-) is implied. The logarithmic 
In (as) part of Eq. (H) has been derived first in 
|17] , |l8|] . The coefficient parameterizes the i- 
loop correction to the Coulomb potential (ai = 
31Ca/9 - 20TFni/9, . . . @). At present, only 
Fade estimates of the three-loop MS coefficient 
03 are available [0. For the bottom and top 
quark case they read 
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The third order corrections to the ground state 
energy proportional to the /3-function only get 
contributions from (i) and (ii) and read 

, 3 
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where — hi(fi/ (CpasiTiq)), C(z) is Riemann's 
^-function and Pi is the {i + l)-loop coeffi- 
cient of the QCD /3 function (/3o = llC^/12 - 
Tpni/3,...). The terms proportional to /Sq in 
Eqs. (^) have been computed first in 

The total result for the third order correc- 
tion to the ground state energy is given by the 
sum of Eqs. (||) and (^) Adopting the choice 
/is = CFOCs(iJLs)'mq one obtains for the bottom 
and top system (for S=l) in numerical form 
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where we have separated the contributions arising 
from as and /3q. The only unknown ingredient in 
our result for SE-i is the three-loop MS coeffi- 
cient as of the corrections to the static potential 
entering Eq. (||). Up to now there are only es- 
timates based on Pade approximation which we 
will use in our analysis. However, our final result 
only changes marginally even for a rather large 
deviations of from its Pade estimate. 

4. PHENOMENOLOGY OF ti THRESH- 
OLD PRODUCTION 

Let us focus on the top quark threshold pro- 
duction, (the bottom quark case is considered 
in p^). The corresponding experimental data 
are not yet available and our goal is to present 
a formula which can be directly used for the top 
quark mass determination from the characteris- 
tics of the cross section of ti production near 
threshold . One has to distinguish the produc- 
tion in e"'"e~ annihilation where the final state 
quark-antiquark pair is produced with S = 1 and 
in (unpolarized) 77 collisions where the dominant 
contribution is given by the S'-wave zero spin final 
state. The nonperturbative effects in the case of 
the top quark are negligible, however, the effect 
of the top quark width has to be taken into ac- 
count properly |^ as the relatively large width 
smears out the Coulomb-like resonances below 
the threshold. The NNLO analysis of the cross 



section shows that only the ground state pole 
results in the well-pronounced resonance. The 
higher poles and continuum, however, affect the 
position of the resonance peak and move it to 
higher energy. As a consequence, the resonance 
energy can be written in the form 
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S^'Er. 
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To compute the shift of the peak position due to 
the nonzero top quark width J'^'Srcs we use the 
result of for the e+e^ ti and 77 ti cross 
sections at NNLO. Numerically for both processes 



5^'E,^^ = 100 ± 10 MeV. 
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This value is rather stable with respect to the 
variation of all input parameters of our analysis. 
To evaluate the perturbative contribution we use 
the result of the previous section up to N^LO. In 
Fig. |l| the perturbative ground state energy in the 
next-to-leading (NLO), NNLO and N^LO is plot- 
ted as a function of the renormalization scale of 
the strong coupling constant for S = 1. One can 
see, that the N'^LO result shows a much weaker 
dependence on /j, than the NNLO one. Moreover 
at the scale /i « 15 GeV which is close to the 
physically motivated scale /Xs, the N'^LO correc- 
tion vanishes and furthermore becomes indepen- 
dent of /i, i.e. the N'^LO curve shows a local 
minimum. This suggests the convergence of the 
series for the ground state energy in the pole mass 
scheme. Collecting all the contributions we ob- 
tain a universal relation between the resonance 
energy in ti threshold production in e+e^ annihi- 
lation or 77 collisions (the result shows very weak 
spin dependence) and the top quark pole mass 
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where the weak nonlinear dependence of Eq. (|6|) 
on mt is taken into account in the second term on 
the right-hand side of Eq. (H). The central value 
is computed for nit — 174.3 GeV. The analysis 
of the theoretical uncertainty in Eq. (H) can be 
found in Due to the very nice behavior of 
the perturbative expansion for the ground state 
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Figure 1. Ground state energy Ef '*" of the ti bound 
state in tlie zero-width approximation as a function 
of the renormalization scale /i for and 5* = 1. 

energy we do not expect large higher order cor- 
rections to our result. As a consequence the top 
quark pole mass can be extracted from the fu- 
ture experimental data on tt threshold produc- 
tion with the theoretical uncertainty of ±80 MeV. 
Thus we find the N'^LO corrections to stabilize 
the behavior of the perturbative series for the 
ground state energy and to be mandatory for a 
high precision theoretical analysis. 
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